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Abstract 

Recently, it has been shown that if we consider the higher derivative correction, 
the viscosity bound conjectured to be 77/s = l/47r is violated and so is the causality. 
In this paper, we consider medium effect and the higher derivative correction simulta- 
neously by adding charge and Gauss-Bonnet terms. We find that the viscosity bound 
violation is not changed by the charge. However, we find that two effects together 
create another instability for large momentum regime. We argue the presence of 
tachyonic modes and show it numerically. The stability of the black brane requires 
the Gauss-Bonnet coupling constant A(= 2a' /t^) to be smaller than 1/24. We draw 
a phase diagram relevant to the instability in charge-coupling space. 



1 Introduction 



After the discovery of consistency of AdS/CFT [1-3] result and that of RHIC experiment 
on the viscosity/entropy-density ratio [4-6], much attention has been drawn to the calcu- 
lational scheme provided by string theory. Some attempt has been made to map the entire 
process of RHIC experiment in terms of the gravity dual [7] . The way to include chemical 
potential in the theory was figured out in [8,9]. Phases of these theories were also discussed 
in D3/D7 setup [10-12] as well as in D4D8D8 [9]. 

More recently, it had been conjectured that the viscosity value of theories with gravity 
dual may give a lower bound for the i]/ s = l/An for all possible liquid [13]. However, the 
authors of [14] and [15] showed that if we consider the stringy correction to a' order, the 
viscosity bound is violated and causality is also [16] violated as a consequence (See also for 
more recent paper [17]). 

The a' terms are also related to the (in) stability issues of black holes. The instability of 
£)-dimensional asymptotically flat Einstein-Gauss-Bonnet black holes has been discussed 
by several authors [18, 19]. Their results show that for the gravitational perturbations 
of Schwarzschild black holes in D — (from 5 to 11) Gauss-Bonnet gravity, the instability 
occurs only ior D — 5 and D — 6 cases at large value of a' [19]. 

In this paper, we add charge together with the Gauss-Bonnet term, and calculate 
rj/s and consider the stability issue including the causality violation. We find that the 
viscosity bound violation is not changed by the charge. However, we find that for large 
momenta regime, there exists a new instability due to the charge effect. The linearized 
perturbation has a negative frequency squared signaling an instability. We draw the phase 
diagram relevant to the instability. The stability of the black brane requires A < 1/24. 
We emphasize that the new instability present only if both charge and Gauss-Bonnet term 
present. 

The rest of the paper goes as follows. In section 2, to set up we give a briefly review on 
the thermodynamic properties of Reissner-Nordstrom-AdS black brane solution in Gauss- 
Bonnet gravity. In section 3, the Gauss-Bonnet correction to rj/s is calculated via Kubo 

formula and its charge dependence is given in an explicit form. In section 4, we study the 
causality violation problem for charged black branes and reproduce the results found in 
Ref. [16]. In section 5, we discuss the stability of Reissner- Nordstrom- AdS black branes in 
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Gauss-Bonnet gravity. Conclusions and discussions are presented in the last section. 



2 Reissner-Nordstrom-AdS black brane in Gauss- 
Bonnet gravity 

The thermodynamics and geometric properties of black objects in Gauss-Bonnet gravity 
were studied in several papers [20-25]. In this section, we mainly review the basic features 
of charged black holes in Gauss-Bonnet gravity. Further details can be found in [23] . 

Wc start by introducing the following action in D dimensions which includes Gauss- 
Bonnet terms and U{1) gauge field: 

/ = j^^Jd''xV^(^R-2A+a' {R^.p^R^""'"' - ^R^.R'"' + i?^) -477(7^ F^.F''^) , (2.1) 

where a' is a (positive) Gauss-Bonnet coupling constant with dimension (length)^ and 
the field strength is defied as F^i,{x) — d^A,^{x) — d,^A^{x). The corresponding Einstein 
equation leads 

R^.u - Iqi^^R + Qt^uA = SttGd [p^pF^.g'"' - ^g^.Fp^F^') + T^^, (2.2) 



where 



T'^ = a' 



-j9vi.v (^RaPpaR"'^'"^ — 4:Ra/3R"'^ + R^^ — '^RRpu + ^RppRy^ 



_|_4 R R P" — OJ^ G P^l 



(2.3) 



The charged black hole solution in D dimensions for this action is described by [23] 

ds^ = -H{r)N^dt^ + H-\r)dr'^ + —hijdx'dx^, (2.4a) 



P 

Q_ 

47r(D - 3)r^ 



with 



^, , , ( / 4a/ mP q^P \\ , (L> -!)(£>- 2) 
^W = ^o + ^^l-yi-^(l-^ + ^jj, A = -^ A 1, 

where a and a' arc connected by a relation a = {D — A){D — 3)a' and the parameter / 
corresponds to AdS radius. The constant A^^ will be fixed later. Note that the constant 



value of ko can be ±1 or and hijdx^dx^ represents the line element of a (D— 2) -dimensional 

hypersurface with constant curvature {D — 2){D — 3)ko and volume Vd-2- The gravitational 

mass M and the charge Q are expressed as 

(D - 2)Vd^2 
M = m, 

27r(D-2)(D-3) 2 
W — ^ Q ■ 

Ltd 

Taking the limit a' — with = 0, the solution may correspond to one for Reissner- 
Nordstrom-AdS (RN-AdS). The hydrodjTiamic analysis in this background has been done 
in [26]. 

In the following, we mainly focus on five-dimensional case with ko = 0. Defining 
A = a/l'^{= 2a'/P), the function H{r) becomes 

^2 



H{r) 



2XP 



[2.5) 



where r+ and r_ correspond to the outer and the inner horizons, respectively, and —r^ = 
t\ + r^. The constant A^^ in the metric f l2.4ap can be fixed at the boundary whose 
geometry would reduce to fiat Minkowski metric conformaly, i.e. ds^ oc — c^dt^ + dx^. On 
the boundary r ^ cxd, we have 

H{r)N^ - ^, 

so that A^^ is found to be 

iV2 



^(l + v/r^). (2.6) 

Note that the boundary speed of light is specified to be unity c = 1. 

We shall give thermodynamic quantities of this background. The temperature at the 
event horizon is defined as 

T= ' (27) 
27r^/^ dr 2ttP \ rl J ' ^ ' ^ 

The black brane approaches extremal when q'^P/r^ 2 (i.e. T ^ 0). The entropy of 
RN-AdS black holes with Gauss-Bonnet terms can be obtained by using 5* = —dF/dT, 
where F is the free energy. After Wick rotation i.e. t ir, the free energy can be obtained 
from the action I in (12.11) through F = —TL By using an identity which is derived by 
taking trace over the equation (12.21) . 

a' {R^ - m^.R'^'' + R^.p^R'^'^P'^) +3R- lOA - AnG^g^'-gP'^F^pF,^ = 0, 
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the action (12. ip reduces to be on-shell, 

1 



n , dr d-'x^l -2R + 8A 



1 V3N 

levrGs TP 
1 V3N 



drr 



-{r'H{r)) -- 



167rG, T/3 



2 r-'H r] 



r+ 
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where we used the scalar curvature i? = — {r^H{r))" /r^ Q. Divergent terms arise in the 
action and we can regulate the result by subtracting the action of the Gauss-Bonnet- 
modified pure AdS space, which is obtained by setting r± = and tq = in H{r){no 
horizons), that is to say 



-^AdS- 



GB 



00 1-13' 



P'VsN 



16nG. P 



2 r-^Hir 



12 



where we assign a temperature P' to AdS space with Gauss-Bonnet terms which is [15,27] 



^BH \ 1/2 

ytt 

„AdS 

ytt 



^BH \ 1/2 

ytt 



,AdS 



Then we find 



AI = 1 -I 



1 UriV-^N V^Nrl 



AdS-GB 



levrGs TP 
Finally, we obtain the entropy densitjC, 

1 dF 1 d{TAI) 1 r\ 



2G.P 



s = — 



dT V:, dT 



4G. P ' 



(2. 



3 Viscosity to entropy density ratio 



Before considering a linearized perturbative theory in this background, we shall summarize 
the basic procedure to calculate Green function and the shear viscosity in Minkowski 
spacetime [28]. We work on the five- dimensional background, 

ds = gmndx"^dx'^ + guudu^, 



*Note that when fco ^ 0, the scalar curvature is R — ~ (r^H{r))" /r^ + 6fco/r^. 

t For more general case, if A; 0, the Bekenstein-Hawking area law is broken and the resulted entropy 



is given by (l + g-f^^o) (see [22] 



where x'" = [t, x, y, z) and u are the four- dimensional and the radial coordinates, respec- 
tively. We refer the boundary as m = and the horizon as u = 1. A solution of the 
linearized equation of motion may be given by, 

where the function fk{u) is normalized such that /a;(0) = 1 at the boundary. An on-shell 
action might be reduced to surface terms in four dimensions by using the equation of 
motion, 

^[0^°^] = / (3.1) 

where the function Q{k,u) can be written in terms of f±k{u) and duf±k{u). Accommo- 
dating Gubser-Klebanov-Polyakov/Witten relation [2,3] to Minkowski spacetime. Son and 
Starinets arrived at the following formulation for the retarded Green function: 

G{k) = 2g{k,u) , (3.2) 

ji=0 

where the incoming boundary condition at the horizon is imposed. In this paper we 
consider the tensor type perturbation in the background. By using an obtained retarded 
Green function, one can estimate the shear viscosity rj via Kubo formula, 

^ _ In.(G(^.0)), ,3.3) 

Now let us proceed to calculate the shear viscosity by using Green function in our 
background. As we see above, it is standard to introduce new dimensionless coordinate 
u = r'^/r'^. The five-dimensional metric with /cq = in fl2.4ap is then deformed into 

where 



f{u) = — 1- v/l-4A(l-M)(l + M-aM2) 
2A L 

and we denote a = q'^P/r'^ and = In this coordinate system, the event horizon of 

the black brane is at m = 1, while m = is the boundary of the AdS space. 

We now study small metric fluctuation hy{t,z,u) = (f){t,z,u) around the black brane 
background of the form 

2 _ -fiu)N'^dt'^ + dx^ + 2(j){t,z,u)dxdy Pdu^ 
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By considering the spin under the 0(2) rotation in (x, |/)-plane, gauge perturbations would 
be decoupled within this tensor type perturbation. Using Fourier decomposition 

we can obtain the following linearized equation of motion for from the equation fl2.2l) : 



= <l>"{u) + J^'P'H 



-m - " - ''t 'fr'""{':»::i'r::^'"'»'^(n). 0.6) 



where 



uN^P{uy^ ' uf{u)[l + 2\u'^{u-^f{u))'] 

g{u) = u-'f{u) [l + 2\u\u-^f{u))'], 
1% - Pb^ 

UJ = OJ, k = K, 

2 ' 2 ' 

and the prime denotes the derivative with respect to u. 

Let us solve the equation of motion (13.61) in hydrodynamic regime i.e. small u and k. 

We first impose a solution as 

(^{u) = {l-urF{u), (3.7) 

where F{u) is a regular function at the horizon m = 1, so that the singularity at the 
horizon might be extracted. Substituting this form into the equation of motion, we can fix 
the parameter u as u = ±ic<j/(47rT) where T is the temperature. We here choose 

UJ 

V = —i , 

47rT' 

as the incoming wave condition. In order to get the shear viscosity via Kubo formula (13. 3p 
by using Green function (13. 2p . it might be sufficient to consider series expansion of the 
solution in terms of frequencies up to the linear order of uj{= iAnTu), 

F{u) = Foiu) + pF^{u) + 0{p\ k^). (3.8) 

The equation of motion (13. 6p becomes the following form up to 0{v), 

[g{u)F'{u)\ - V (y^^7h) F{u) - -^g{u)F'{u) = 0. (3.9) 

Substituting the series expansion (13.80 into the equation (13.90 . one can get the equations 
of motion for Fq{u) and Fi{u) recursively. From Olu^) in the equation (13.91) . the equation 
for Fq{u) is obtained as 

[giu)F;,iu)]' = 0, (3.10) 
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and can be solved as 

where Ci is an integration constant. Regularity of -Fo('^^) at the horizon implies that Ci 
must be zero as g{u) goes to zero at the horizon. Therefore, Fo{u) is a constant, 

Fo{u) = C, (const.). (3.11) 

The solution for Fi{u) can be obtained from equation (13. 9p at O^u^), 

[9{u)F[{u)]' - (^-^g{u)^' = 0. (3.12) 
Integrating the above equation we get 

F[{u) = -^ + -^ (3.13) 

1-u g[u) 

The integration constant C2 can be fixed by the regularity condition of Fi{u) at the horizon. 
At the horizon m = 1, the function g[u) behaves as 

g{u) = ((2 - a)(l - 2A(2 - a))) (1 - n) + 0((1 - uf). 

Then, the regularity condition at m = 1 implies 

C2 = -(2-a)(l-2A(2-a))C. (3.14) 

The remaining constant C is estimated in terms of boundary value of the field, 

\im (j){u) = 

SO that we could fix 

C = 0W(^l + C(z/)). (3.15) 

Now we shall calculate the retarded Green function. Using the equation of motion, the 
action reduces to the surface terms. The relevant part is given as 



n=l 



u=0 



(3.16) 



Near the boundary u = e, using the obtained perturbative solution for 0(m), we can get 

_ _^(2-a)(l-2A(2-a)) ^ 
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Therefore we can read off the correlation function from the relation (13.21) . 

Gxy xy{u, k) = -^^Ye^ ) " ^^^^ " + ^^-^^^ 

where we subtracted contact terms. Then finally, we can obtain the shear viscosity by 
using Kubo formula (13.31) . 

The ration of the shear viscosity to the entropy density is concluded as 

One can see explicitly that the conjectured viscosity bound can be violated for some value 
of A and a. Figure [1] demonstrates that for fixed value of r]/s, as the couping constant A 
increases, a also increases. The shear viscosity approaches zero as (A, a) (1/4,0) and A 
is thus bounded by 1/4. When a = (no charges), i]/s = (1 — 4A)/(47r), we recover the 
result in Ref. [15]. It is also worth noting that for extremal case (a = 2), the ratio of the 
shear viscosity to entropy density receives no corrections from Gauss-Bonnet terms. 




0.00 0.05 0.10 0.15 n.2n o.is 



Figure 1: Shear viscosity to entropy density ratio as a function of a and A. The lines correspond 
to rj/s = 0.07, 0.06, 0.01, respectively, from top to bottom. 
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4 Causality violation 



It was shown that the causahty could be violated if one introduced Gauss-Bonnet terms [15, 
16]. We here analyze an effect of the charge to this issue. 

Due to higher derivative terms in the gravity action, the equation (13.61) for the prop- 
agation of a transverse graviton differs from that of a minimally coupled massless scalar 
field propagating in the same background geometry. Writing the wave function as 

^{x,u) = 6-'^'+''''+'''-'', (4.1) 

and taking large momenta limit oo, one can find that the equation of motion (13. 6p 

reduces to 

n'^gf, ^ 0, (4.2) 

where the effective metric is given by 

Note that can be interpreted as the local speed of graviton: 



N^fiu) 1 - 2Au2(2n(M-i/(u))" + 3(u-i/(n))' 



1 + 2Am2(m-i/(m))' 
We can expand near the boundary u = 0, 

5,. . 2(1 + a) 1 



(4.4) 



cj-l= --(l + a)+ ^ / u' + 0(u'). (4.5) 

As we will see below, the local speed of graviton should be smaller than 1 (the local speed 
of light of boundary CFT). We require 

- - + < 0. (4.6) 

2 1-4A 2v/r^4A- ^ ^ 

The above equation leads to A < 0.09 without any charge dependence. This is the same 

result with neutral black holes in Gauss-Bonnet theory [16]. For charged black branes, the 

speed of graviton is also smaller than 1 in the range A < 0.09. 

Now let us study in the regime A > 0.09, to see how the causality is violated in the 

boundary theory. We follow the discussion in the papers [15,16]. By using the identification 

— — = g'^^^'^kjj, we may rewrite the equation (14. 2 h into one for geodesic motion, 
ds 

dx'^ dx^ 
ds ds^^" 



gfu = 0. (4.7) 
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If we consider uj and q as conserved integrals of motion along the geodesic with 



from (14. 2p and (14.71) . we can obtain 



(4.9) 



Ffo ds / P f 

We can simplify the above equation by recalling s as s = ks/N and noting that u = 

(|)^a.^-4 a^^^. (4.10) 

The geodesic equation determines the radial motion of a test particle with energy in an 
effective potential c^. One can now infer that, from Figure [21 geodesic line which starts 
at spatial infinity (the boundary) can bounce back to the boundary. In other words, the 
turning point appears at 

= c'Juo). (4.11) 




Figure 2: The hump (I) signifies the causality violation, while the well (II) indicates that the 
black brane is unstable. 

For a light-like geodesic line starting from the boundary and bounced back at the boundary, 
we find 

Az = 2 ^^dr = 4-r / dr ' (4.12b) 
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As pointed out in Ref. [16], one may find microcausality violation in the boundary CFT 
when the a bouncing graviton geodesic satisfies > 1. From fl4.12ap and fl4.12bp . we can 
see that if we tune Cg{uo) to be c^^max, we then have 

Az 

Cg,m^^ > 1. (4.13) 

Since near the boundary Cg can be greater than 1, the propagation of signals in the boundary 
Az 

theory with speed might become superluminal. Actually, metastable states can live 

Az 

in the well between cl^.^^ and the boundary and — — corresponds to the group velocity 

^' At 

of metastable quasiparticles. According to the standard procedure to analyze quasinormal 
modes, we rewrite the wave function in a Schrodinger form, 

where ^/'(r(r*)) and the potential is defined by 



^ = Kir)^, Kir) ^ J = 1 - Xbr^^'^'^ ^^'^^ 



u ^f{u) dr 



V = k'^cl + Vi{r), Vi{r) = NH^b'^ 



d\nK{r)\ d f dhiK{r 



fir) ^ +/(r)^ /(r 



dr / dr \ dr 



Following the procedure of Ref. [16], one can find that the group velocity of the graviton 
is given by 

Therefore, signals in the boundary theory propagate outside of the light cone. Again, we 
confirm that microcausality violation happens in the CFT. One may expect that when 
A < 0.09, the theory with Gauss-Bonnet corrections is safe and consistent. 

5 Instability 

Apart from the causality violation, for RN-AdS black brane in Gauss-Bonnet theory, the 
charges give the instability of the black brane within the window of < A < 0.09. 

From Figure [HI we can see that the Schrodinger potential develops a negative gap near 
the horizon. We will now show that in the large momentum limit, the negative- valued 
potential leads to instability of the black brane. In the large momenta limit k^^ — oo, the 
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Figure 3: Schrodinger potential V{r) as a function of r for A = 0.2, a = 1.7 and k = 500. Near the 
horizon, the potential develops a negative- valued well and then negative-energy bound states will 
appear. Near the boundary, the potential develops a hump which corresponds to superluminal 
propagation of metastable quasiparticles. 

dominant contribution to the potential is given by k'^c^- For near extremal cases, can be 
negative near the horizon and V ~ k'^c^ can be deep enough (see Figure [3]). Thus bound 
states can live in the negative-valued well. The negative energy bound state corresponds 
to modes of tachyonic mass on Minkowski slices [29] and signals an instability of the black 
brane [18, 19]. Let us expand in series of (1 — m), 

(2 - a) f 1 + 4A - 14aA - 32X'^ + 32aA2 - Sa^xA 
4- ^ iT^IxT^) + (5.1) 

Since < a < 2, and < n < 1, will be negative, if 

(1 + 4A - 14aA - 32A2 + 32aX'^ - Sa^A^) 
(1 - 4A + 2aA) 

From the above formula, we find the critical value of A, 



< 0. (5.2) 



2 - 7a + V3V12 - 20a + 19a^ 
= 8(^32? • ^'-'^ 

Above the line of Ac, can be negative (see figure 2). The minimal value of Ac can be 
obtained in the limit a — 2, 

Ac, mill = (5.4) 

Figure H] shows us that the two lines Ac (a) and A = 0.09 separates the physics into 
four regions in (a. A) space. The physics in region I so far can be consistent. In region 
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Figure 4: Phase diagram for the instabihty in a-X space. Region I: There are no causahty 
violation and bulk bound states. Region II: Causality violation can happen, but bound states do 
not appear in the bulk. Region III: Both causality violation and instability happen in this region. 
Region IV: There is no causality violation, but the black brane is unstable. 

II, causality violation can be found. In region III, causality violation as well as unstable 
quasinormal modes (QNMs) appear. In region IV, we can only find unstable QNMs. 
Figures [5] and [6] show us explicitly the behaviors of in different regions. In order to 
demonstrate that the peculiar feature of < signals the instability of the RN-AdS black 
brane in Gauss-Bonnet theory, we solve the Schrodinger equation f l4.14p with negative- 
valued potential numerically and find some unstable QNMs (see Table 1). From Table 1, 



Table 1: Unstable QNMs for charged GB black brane perturbation of tensor type. 



A 


a = 1.9 


a = 1.7 


a = 1.4 


a = 1.1 


a = 0.8 


a = 0.5 


0.2 


258.2823i 


255.9158i 


245.6454i 


223.4933i 


179.9353? 


94.4075? 


0.15 


158.0242? 


154.8652i 


141.6929i 


113.4577i 


60.5891? 




0.1 


79.3897i 


75.4296i 


59.4755? 


26.9101Z 






0.07 


33.9841Z 


29.3078i 










0.05 


6.6761i 













we can find that the real part of uj is vanishing, while the imaginary part of oo is positive. 
Inserting the values of u into the equation (14. ip . one can see that gravitational instability 
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Figure 5: as a function of u. Line I 
describes the behavior of in region I for 
A = 0.05, a = 0.2. Line II shows that > 1 
at some value of u near the boundary for 
A = 0.14, a = 0.4. 



Figure 6: as a function of u. Line III: has 
a hump greater than 1 and a negative-valued 
well (We set A = 0.15, a = 1.7). Line IV: 
only has a negative-valued well (A = 0.089, 
a = 1.99) 



grows as time goes on and then the black brane becomes unstable against gravitational 
perturbation. The numerical analysis also indicate that the black brane becomes stable 
under gravitational perturbation, when we restrict A to be A < 1/24. 

6 Conclusions and discussions 

In summary, we have computed the charge dependence oirj/ s for Gauss-Bonnet theory. We 
have taken RN-AdS black brane solution into Gauss-Bonnet gravity and used the Kubo 
formula to compute the viscosity of the dual boundary theory. The ratio of the shear 
viscosity to entropy density was found to he r]/ s = {l — 4A(1 — a/2) j / (47r), which violated 
the conjectured viscosity bound for non-extremal RN-AdS black brane in Gauss-Bonnet 
gravity. However, for extremal case, the conjectured lower viscosity-entropy density bound 
l/47r can be recovered. 

The causality violation and the instability of charged black brane were also analyzed 
in this paper. We have confirmed the results found in previous work that when A > 0.09, 
causality violation happens in the boundary CFT [16]. It is interesting to notice that 
charges introduce instability of the black brane even in the range < A < 0.09. Therefore, 
to avoid causality violation and instability for any charge, we suggest to restrict the value 
of A to be A < 1/24 ~ 0.04. 
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Our final remark is about the interpretation of tlie cfiarge effect. If one introduces tfie 
bulk filling D-branes, one can consider the two kind of Maxwell fields, one for i?-charge 
U{1) and the other for the baryon charge for the brane gauge field [30]. Since both fields 
couple to the gravity in the same way, one can consider the charge in our analysis either 
as the i?-charge or baryon charge. For the latter case, we can interpret the charge effect 
as the effect of finite baryon density. 
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